Density Wave States of Non-Zero Angular Momentum 
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We study the properties of states in which particle-hole pairs of non-zero angular momentum 
condense. These states generalize charge- and spin-density-wave states, in which s-wave particle- 
hole pairs condense. We show that the p-wave spin-singlet state of this type has Peierls ordering, 
while the d-wave spin-singlet state is the staggered flux state. We discuss model Hamiltonians which 
favor p-wave and d-wave density wave order. There are analogous orderings for pure spin models, 
which generalize spin-Peierls order. The spin-triplet density wave states are accompanied by spin- 
1 Goldstone bosons, but these excitations do not contribute to the spin-spin correlation function. 
Hence, they must be detected with NQR or Raman scattering experiments. Depending on the 
geometry and topology of the Fermi surface, these states may admit gapless fermionic excitations. 
As the Fermi surface geometry is changed, these excitations disappear at a transition which is third- 
order in mean-field theory. The singlet d-wave and triplet p-wave density wave states are separated 
from the corresponding superconducting states by zero-temperature C'(4)-symmetric critical points. 



PACS numbers:71.10.Hf, 75.10.Lp, 75.30.Fv, 71.27.-(-a, 



I. INTRODUCTION 

In recent years, a number of materials have been un- 
covered in which the competition between an effective at- 
tractive interaction and short-range repulsion appears to 
lead to the formation of superconducting states in which 
the Cooper pairs have non-zero relative angular momen- 
tum. In this paper, we suggest that such competition can 
also lead to density-wave states formed by the conden- 
sation of particle-hole pairs of non-zero relative angular 
momentum. These states generalize the familiar charge- 
and spin-density-wave states, in which s-wave particle- 
hole pairs condense. We discuss several different possible 
ordering schemes, the types of interactions which favor 
them, their physical properties, and their possible rele- 
vance to experiments. 

Several such states are already commonly known by 
other names, as we will show below. The singlet / = 1 
density-wave state is simply the Peierls state (or bond- 
ordered wave) , while the singlet I = 2 density- wave state 
is known as the staggered flux state of [Q-^ . However, the 
triplet analogues of these states have not been discussed. 
Since the triplet analogues of these states break spin- 
rotational invariance, they have S = 1 Goldstone boson 
excitations. However, the ground state does not have a 
non-zero expectation value for the spin at any wavevec- 
tor. Hence, as we will see, these Goldstone bosons cannot 
be detected in experiments which couple simply to the 
spin density, such as neutron scattering or NMR. Instead, 
Raman scattering or NQR are necessary to couple to the 
Goldstone bosons of these more subtle types of ordering. 
More generaly, s-wave probes cannot couple directly to 
the orders discussed here; instead, local probes or those 
which couple to higher powers of the order parameter are 
necessary. 

The p-wave and d-wave density wave states are favored 
by the same types of interactions which favor the s-wave 



state - i.e. the CDW. However, they evade interactions 
which disfavor CDW order. Similarly, they are favored 
by superconductivity-favoring pair-hopping terms 
while evading interactions which disfavor superconduc- 
tivity. Hence, they are rather natural candidates for sys- 
tems with competing repulsive interactions. 

As in the case of higher angular momentum supercon- 
ducting states, there is the possibility of gapless excita- 
tions since the order parameter can have nodes on the 
Fermi surface. To consider one consequence of this, sup- 
pose that the shape of the Fermi surface is such that the 
nodes of the order parameter do not lie on the Fermi 
surface. Let us distort the shape of the Fermi surface by, 
say, changing the anisotropy between the hopping param- 
eters, which can be done by applying uniaxial pressure. 
At the mean-field level, a third-order phase transition 
can occur at which gapless excitations appear. After this 
point, the system remains critical as a result of the node. 

The analogy with supercondutivity can be taken a step 
further by combining density wave order with supercon- 
ducting order in a pseudospin SU (2) triplet following 
Yang [|j and Zhang At a critical point between the 
two types of order, this pseudospin SU{2) could become 
exact, giving - together with SU{2) spin symmetry - an 
0(4)-invariant critical point. We discuss the possible rel- 
evance of such a critical point to the pseudogap regime 
of the cuprate superconductors. 

Particle-hole condensates with non-zero angular mo- 
mentum were considered in the context of excitonic in- 
sulators by Halperin and Rice [Q. They were redis- 
covered in the context of the mean-field instabilities of 
of extended Hubbard models by Schulz ||l^ and Ners- 
esyan [|l2|,|l3| and collaborators. At around the same 
time, Kotliar [Q| and Marston and Affleck ||^ found the 
staggered flux state as a mean-field solution of the Hub- 
bard model. However, it was apparently not recognized 
that the singlet d,j.2 _y2 density- wave state is the same as 
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the staggered flux state. More recently, this state 14 -T^] 
and a related variant have been discussed in the 

context of the cuprate superconductors. A version of this 
state (see the comments in the concluding section) has 
appeared in mean field analyses of an SU (2) mean field 
theory of the i - J model The Nodal Liquid 



state of [21 23 also bears a family resemblance to the 



staggered flux state; we will return to the relationship 
between these states in the concluding section. 



II. ORDER PARAMETERS AND BROKEN 
SYMMETRIES 

We define the different possible density-wave orderings 
by analogy with the more familiar superconducting case. 
Consider a system of electrons on a square lattice of side 
a. A superconductor is defined by a non- vanishing ex- 
pectation value of 



(1) 



A triplet superconductor is characterized by the expec- 
tation value 



{ipaik,t)'ipf3{-k,t)) = A{p) ■ a^e^p 



(2) 



Fermi statistics requires that A(p) be odd in p. p-wave 
superconductors can have the components of A(p) chosen 
from sinfcj^a, sinfcj,a, or sin fc^a i * sinfcj^a. For instance, 
a Px superconductor with all spins polarized along the 3- 
direction will have Ai -\-i/S.2 ^ and A3 = Ai — iA2 — 0: 

{■ijja{k,t)'ijjf3{-k,t)) = Ao (sinfc^jfl) ct+ ''e^p (3) 

Spin-polarized Py and Px + ipy superconductors have 
smkxa replaced, repectively, by sinfcj^a and svukxa + 
isinkyQ. The analog of the A' phase of ^He has equal 
numbers of tt a-nd |J, pairs: 

Aq (sin k^a cr\^ + sin ky a a'^ e^fj (4) 

An unpolarized superconductor of ti pairs has A3 ^ 
and Ai = A2 = 0: 

{'ipa{k,t)'ipf3{-k,t)) = Ao (sinfc^a) al '^e^p (5) 

As in the polarized case, unpolarized Py and Px + ipy 
superconductors have sinfe^^a replaced, repectively, by 
sinfcya and sinfca;a-|- jsinfcyO. In principle, more compli- 
cated order parameters are possible, with all components 
of A taking non- vanishing values. If any component of 
A(p) is not real, time-reversal symmetry (T) is broken. 

A d-wave superconductor must be a spin-singlet super- 
conductor. A dx2-y2 superconductor has 

{4>a{kjt) 'ipf3{—k,t)) — Aq (cos kxU — cos kyu) (6) 



while a dxy superconductor has cos fc^^a — cos kya replaced 
bysinfcj^a sin/cyO. A (i2,2_y2 -l-irf^j, superconductor breaks 
T with the order parameter: 

{',p^{k,t)^p{-k,t)) = 

Aq (cos kxtt — COS kya + i sin kxO. sin kyo) eap (7) 

We can define analogous orders for density-wave states. 
However, the spin structures will no longer be determined 
by Fermi statistics. Let us first consider the singlet order- 
ings. A singlet s-wave density wave is simply a charge- 
density-waveQ: 

{rHk + Q,t)4^0{k,t))^<^>QS'^ (8) 

A singlet px density-wave state has ordering 

(V-^^Cfc + Q, t) 4>p{k, t)) = $Q sin kxa S'^p (9) 

The singlet px -f ipy density-wave states are defined by: 

(■(/'"■'■(A: + V/3(fc,t)) = $Q {sinkxa + i sin ky a) 5^ 

(10) 

Similarly, the singlet dx2^y2 density- wave states have 

{tp"'^{k + Q,t) iif3{k, t)) = <I>Q (cosfcajfl - coskyo) 8'p 

(11) 

while the singlet dx^^y^ -\- idxy density- wave states have 

{r\k + Q,t)Mk,t))^ 

<f>Q (cos kxtt — COS kya + isinkxa sinkytt) 6^ (12) 

These states belong to a class of states of the form: 

(i^^Hk + Q, t) Mk, t)) = $Q /(fc) 5^ (13) 

f{k) is an element of some representation of the space 
group of the vector Q in the square lattice. In this paper, 
we will focus primarily on the cases f{k) = sinkxa and 
f{k) = cos kxa — cos kyQ, but f{k) could be an element 
of some larger representation. The s-wave (or extended 
s-wave) cases, f(k) = |/(fc)|, are the usual charge-density 
wave states. 

Q is the wavevector at which the density-wave or- 
dering takes place. It may be commensurate or in- 
commensurateQ. For commensurate ordering such that 
2Q is a reciprocal lattice vector, e.g. Q ~ (7r/a,0) or 



* Extended s-wave is also possible. 

^In this paper, we will take commensurate to mean the 
situation in which 2Q is a reciprocal lattice vector. The 
term 'incommensurate' will actually include higher-order 
commensurability. 
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Q = (jr/ajiT/a), we can take the hermitian conjugate of 
the order parameter: 

(V/^(fc,i)^„(fc + g,i)) =^*Qnk)s^ 

{^^\k + g + Q, t) ^„(fc + Q, t)) ^ r{k) 51 

'^Qfik + Q)61 = %rik)61 (14) 



Therefore, for Q commensurate 
f{k + Q) 



rik) 



Q 



(15) 



Hence, if f{k + Q) — —f*{k), $q must be imaginary. For 
singlet Px ordering, this wiU be the case if Q = (tt/a, 0) 
or Q = (7r/a, 7r/a). For singlet dx2-y2 ordering, this will 
be the case if Q = (tt/o, n/a). If f{k + Q) = f*{k), 
must be real. For singlet ordering, this will be the 
case if g = (0,7r/a). For singlet dxy ordering, this will 
be the case if g = (vr/a, ir/a). 

For incommensurate ordering, <i>Q can have arbitrary 
phase: the phase of is the Goldstone boson of bro- 
ken translational invariance, i.e. the sliding density-wave 
mode. Impurities will pin this mode - at second order in 
the impurity potential, as in the case of a spin-density- 
wave - so we will not consider it further. 

All of these states break translational and rotational 
invariance. To further analyze the symmetries of these 
states, it is instructive to write these orderings in real 
space. The singlet Px density- waves have non- vanishing 
expectation value: 



(^^"(x, t) ^pix + ax, t) — ^'^"'{x, t) tI)[3{x — ax, i)) 



Q 



JQ-S 



-iQ-x 



51 (16) 



We have only written the modulated term; the . . . refers 
to the uniform contribution coming from the Fourier 
transform of ip'' {k)ip{k). 

Let us consider the commensurate and incommensu- 
rate cases separately. The incommensurate singlet px 
density-wave states completely break the translational 
and rotational symmetries. If = ~^*_q, T is pre- 
served; otherwise, it is broken. The singlet Px + ipy 
density-wave states always break T. The commensurate 
states, on the other hand, break translation by one lattice 
spacing; translation by two lattice spacings is preserved. 
From (|l|), a commensurate singlet px density- wave state 
with g = (7r/a,0) must have imaginary $q: 



(^ip''°'{x, t) ■0/3(5' + ax,t) — ^'^"{x, t) ijjfj{x — ax,t)) 



"5$ (17) 



The singlet state of this type breaks no other symmetries; 
it is usually called the Peierls state or bond order wave. 
If g = (0,7r/a), $Q must be real. 

(ip''°'{x, t) ipfi^x + ax, t) — ip'^°'{x, t) ipp^x — ax, t)) = 



...-i\^Q\e'^-'51 (18) 

As a result of the i, the Q — (0, n/a) singlet px density- 
wave states break T . However, the combination of T and 
translation by an odd number of lattice spacings remains 
unbroken. The same is true of the commensurate singlet 
Px + ipy density- wave states. Examples of commensurate 
and incommensurate singlet Px and Px + ipy density- wave 
states are depicted in figure |l|. 



(a) 



(b) 



(c) 



(d) 



FIG. 1. (a) Q — (7r/a,0) px density-wave state. (b) 
Q — {OjTv/a) Px density-wave state, (c) Q = {n/a,0) Px +ipy 
density-wave state. (d) Incommensurate Px density-wave 
state. Arrowless lines are bonds where the kinetic energy 
is large but there is no net current. Line thickness indicates 
bond strength. Arrowed lines denote currents. 

The singlet dx^^y^ density- wave states have non- 
vanishing expectation value: 

{^p^" (x, t)ij[3{x + ax,t) + ip^" {x, t) ippix - ax,t)) 
- (^^"(x, t) ^0{x + ay, t) + , t) ^0{x - ay, t)) = 

■ • • + ^ (*Q e''^ " + (19) 



The incommensurate singlet dx^-y^ density- wave states 
will preserve T if <I>q = $1q; otherwise, they break 
T. The same is true of the incommensurate singlet dxy 
density-wave states: 

{ip^"{x, t) ippix + ax + ay, t)) + 
(■0 '''"(:?, t) ^p{x-ax - ay,t)) - 
(ip^"{x, t) i'pix — ax + ay, t)) — 
(^tp''°'{x,t)^f}{x + ax ~ ay,t)') — 

f $_Q e-''^-^) 61 (20) 



Incommensurate singlet dx-^-y-^ +idxy density- wave states 
necessarily break T: 
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{tp'^°'{x, t) tp/six + ax, t) + ip''°'{x, t) ippix - ax, t)) 
- ( V-t" (x, t) ^0 {x + ay,t)+ xjj^" {x, t) i^p {x-ay,t)) 

+i (^p'^"'{x, t) ippix + ax + ay, i)) 

+i t) ippix — ax — ay, i)) 

—i {^p'^°'{x, t) ■0/3(5' — ax + ay, tfj 
—i (tP^°'{x, t) ippix + ax ~ ay, t)) = 

■■■+{l-l){^Q + '^-Q ^'"^'') (21) 

The commensurate Q = {'r:/a,'K/a) singlet d^^-y^ 
density- wave states must have imaginary $q: 

(^^■"(f, t) %l:p{x + ax,t)) + (0^"(x, t) iPfiix- ax,t)) - 
{^^''ix,t) Mx + ciy,t)) + {^^'^ix,t) i;p{x-ay,t)) = 

■■■ + \\'^Q\e''^-'S"p (22) 

As a resuh of the i, the singlet d^^-y^ density- wave 
breaks T as well as translational and rotational invari- 
ance. The combination of time-reversal and a translation 
by one lattice spacing is preserved by this ordering. The 
commensurate Q = (7r/a,7r/a) singlet d^^-y^ density- 
wave state is often called the staggered flux state. There 
is also a contribution to this correlation function coming 
from -ip^ {k)ip{k) which is uniform in space (the . . .); as 
a result, the phase of the above bond correlation func- 
tion - and, therefore, the flux through each plaquette - 
is alternating. The commensurate Q = (7r/a,7r/a) sin- 
glet dxy must have real ^q; therefore, it does not break 
T. On the other hand, the singlet d^i^y-i + idxy state 
does break T. Note that the nodeless commensurate 
singlet dx2-y2 + idxy density- wave state does not break 
more symmetries than the commensurate singlet dx^-y^ 
density-wave state, in contrast to the superconducting 
case. Examples of singlet dx^-y^, dxy, and dx^-y^ +idxy 
density- wave states are depicted in figure ^. 

We now consider the triplet density-wave states. 
Triplet states all break spin-rotational invariance and, 
therefore, have Goldstone boson excitations. We will dis- 
cuss the experimental consequences of these Goldstone 
bosons later. The triplet s-wave density wave state is 
simply a spin-density-wave. The triplet p-wave and d- 
wave states are characterized by: 

{^''\k + Q,t)Mk,t))^^Q{k)-a'l (23) 

with the components of $Q(fc) chosen from, respectively, 
sin kxtt, sin kya, sin kxa ± i sin kya; and cos kxa — cos kya, 
sin kxa sin kya, cos kxa — cos kya ± i sin kxa sin kya. 

A state in which the particle-hole pairs are polarized, 
which is the most direct analogue of a spin-density-wave 
has ^q(p) of the form <i>^ 7^ 0, <i>Q = $q = 0: 

(0"t(^ + t) ^^(k, t)) = -^Q f{k) al " (24) 



(a) (b) 




FIG. 2. (a) Commensurate d^2_y2 density-wave state, 
(b) Commensurate d^y density-wave state, (c) Commensu- 
rate d^2_y2 + idxy density- wave state, (d) Incommensurate, 
T-preserving d^2_y2. Arrowless lines are bonds where the ki- 
netic energy is large but there is no net current. Line thickness 
indicates bond strength. Arrowed lines denote currents. 

where f{k) is chosen from the above set. Alternatively, 
the particle-hole pairs can be unpolarized, e.g. 

{■il:°''i{k + Q,t)ip,3{k,t)) = (sin kxaal " + i sin kyaaj") 

(25) 

As in the superconducting case, more complicated order 
parameters are possible, with all components of A taking 
non- vanishing values. 

For commensurate ordering, we can follow the same 
logic as in (p^. The phases of the components of $q(p) 
are constrained in the same way as the singlet order pa- 
rameters, as illustrated by the i in front of the second 
term in (|2|). 

The orders discussed here can be generalized to other 
2D lattices and to 3D lattices. The orbital wavefunctions 
sinkxa, etc. will be replaced by representations of the 
point groups of these other lattices. 

To each T-preserving singlet ordering, we can associate 
an ordering of a pure spin model, in the same way that 
spin-Peierls ordering is related to Peierls ordering: 

{iP"\k + Q,t)iP^{k,t)) ^ (^Sik + Q,t) ■ Sik,t)'^ (26) 

These spin orderings are states in which the exchange 
energies are large along preferred directions. These pre- 
ferred diractions oscillate from one lattice point to the 
next with spatial frequency Q. The simplest case is spin- 
Peierls ordering, in which the spins form dimers. Another 
example is dxy ordering of a spin model, which takes the 
form: 
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(^S{x, t) S{x + ax + ay, t)^ + (^S{x, t) S{x — ax — ay, t) 
— (^S{x, t) S{x — ax + ay, t)^ — (^S{x, t) S{x + ax — ay, t) 

in analogy with (^). 

III. MODEL HAMILTONIANS 

We are primarily concerned in this paper with the uni- 
versal properties of the states introduced above. We will 
not attempt to show that particular realistic models of in- 
teracting electrons have p- or d-wave density- wave ground 
states. Rather, we will content ourselves with discussing 
the types of interactions which favor such orders and 
showing that they lead to energetically favorable trial 
variational wavefunctions for some idealized Hamiltoni- 
ans. 

The analog of the BCS reduced Hamiltonian for singlet 
density- wave order is: 



H = 



ij''Hk + Q)Mk)^^Hk')^p{k' + Q) 



(27) 



In the triplet case, we replace the four-fermion operator 
of (0) by 



(28) 



We now introduce the variational wavefunction 



= n + "fe,oV'"^(fc + Q)) 



k.a 



(29) 



Its energy can be minimized if we take 

5*Q/(fc) 



(6(fc)-e(fc + Q))VV|$Q|2(/(fc))^ 



(30) 



in the singlet case and 



a^Qfik) 



(e{k)-eik + Q)f+Ag^\<^Q\\f{k)f 



(31) 

in the triplet case, and require <I>q to satisfy the gap 
equation: 



_fk_ 



ifik)f 



(6(fc)-e(fc + g))W|cOQ|2(/(fc))2 



(32) 



The reduced Hamiltonian has long-ranged interactions, 
so the variational wavefunction is essentially correct. We 
will now show that short-ranged Hamiltonians will in- 
clude terms of the form (|2^), and that the trial wave- 
function (^9|) is reasonable for these short-ranged Hamil- 
tonians. 

Consider, then, the following lattice model of interact- 
ing electrons: 

<ij> i 

+ V niUj 

<ij> 
<i,j> ,<i' ,j> 



(4 



(33) 



The first two terms are the usual hopping, t, and on-site 
repulsion, U of the Hubbard model. The third term is 
a nearest-neighbor repulsion, V . The third and fourth 
terms lead to the correlated motion of pairs of electrons, 
tci hops an electron from i' to j when j is vacated by an 
electron hopping to i. tc2 hops nearest-neighbor pairs in 
the same direction. Terms of this general form have been 
discussed in as a mechanism for superconductivity. 
As we will see below, they not only favor superconduc- 
tivity, but p- and d-wave density wave order as well. 

Fourier transforming the interaction terms into mo- 
mentum space, we see that terms of the form of the re- 
duced interaction (p7|) are, indeed, present: 



■ dPk 




Uij^\ki)i;^{k2)^^Hk^)i;^{ki) 
cos (fcg — k1)(ij 



-kl)a' 

- 2tci {cos {kl -kl)a + cos {k\ - fcf ) a 
-t- 2 cos k^a cos fcf 



fcf a -I- 2 cos kfacos k\a 



2tc2 ( sin kfa sin k^a + sin kfa sin k^c 



X ^"^(fci)Vo(fe2)V'^Ufc3)V/3(M 



(34) 
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Let us now consider various candidate orderings and 
the terms which favor or penahze them. Antiferromag- 
netic order is favored by U but penaUzed by V. Charge- 
density-wave order is favored by V but penahzed by U. 
p and d-wave superconductivity are favored by ic2 and 
tci respectively and penahzed by V. p and d-density- 
wave order are favored by tc2 and td, respectively, and 
are both favored by V. The density-wave states can be 
favored over the others by taking V large. The p- or 
d-wave states can be favored by taking tc2 or td large. 
To be more precise, the mean-field equations for various 
ordered states read: 



A 



if{k)f 



(2^^)' ^ie{k)-e{k + Q)f + AX^\^Q\'ifik)f 



1 



(35) 



where 



XdDW = 8V + 96td 
XpDW ^AV + 16tci + 16tc2 
\cDW = 16^ + 96tei -f 16tc2 
Xaf = 2C/ 



2U 



(36) 



Hence, the singlet d^^^ 
ground state if 



density-wave state will be the 



8ici <U -AV < 

8tc2 <2V + AOtd 



(37) 



while the singlet Px density-wave state will be the ground 
state if 



6V + AOtd <U <2V 
2V + 



f 8ici + 8ic2 



(38) 



assuming that the van Hove singularities are at the antin- 
odes of the order parameters. Otherwise, the p- and d- 
wave density wave states will be favored over somewhat 
smaller regions of parameter space. By including spin- 
dependent interactions such as J Si ■ Si, we can favor the 
triplet p- or d-wave states. Hence, it appears that the or- 
derings discussed in this paper are viable. The detailed 
energetics at large coupling strengths - which surely hold 
in physically interesting systems - are beyond the scope 
of this paper. 



neutron scattering. When another symmetry - in addi- 
tion to translational invariance - is broken, this is easier. 

Let's first consider broken time-reversal symmetry. 
The commensurate singlet dx2-y2 density wave state - 
or staggered flux state [^-^ - breaks T; there is an alter- 
nating pattern of currents circulating about each plaque- 
tte of the lattice. These currents produce an alternating 
magnetic field measurable by /iSR and, in principle, by 
neutron scattering The magnitude of the current 
along a link of the lattice will be: 



10 '"^Ampere x $q 



(39) 



Now, $Q is related to the maximum of the gap ac- 
cording to Ao = g^Q where g is the appropriate cou- 
pling constant. Let us suppose, for the purposes of 
illustration, that the formation of the ordered state 
is driven by XdDW- Then, for XdDW small, $q 
(VArfDw)e"'^™"*'*-^*/^''°"'- Alternatively, we may take 
the high-Tc context as a guideline: observed gaps are 
~ 100 — 300K, while interactions such are ~ leV. In 
this case, we expect 4>q ^ 10^^. This translates to a 
magnetic field at the center of each plaquette on the or- 
der of lOG. The muons in a fxSR experiment might see 
a lowwer field if they sit at points of high symmetry or 
away from the plane. The orbital magnetic moments are 
likely to be dwarfed by local spin moments Incom- 
mensurate ordering may or may not break T; if it does, 
the above analysis applies. 

In the <I>^ ^ 0, <i>Q = (f>Q = triplet dx2-y2 density 
wave state, there are counter-circulating currents of up- 
and down-spin electrons. These currents cancel, so there 
is no net current circulating about each plaquette, but 
there is an alternating pattern of spin currents circulating 
about each plaquette. The checkerboard pattern of spin 
currents will generate, via the spin-orbit coupling, 



Hso 



d^k 



2 f^^(l) 



(27r)2 (27r) 



2k + qj 
i;^^{k + q)a/Mk) 



(40) 



a quadrupolar electric field which is, in principle, mea- 
surable in NQR experiments. With the above estimate 
of the current, a nucleus with a non-zero quadrupole mo- 
ment would have an induced splitting of order lOHz. 

We now turn to broken spin-rotational invariance, 
characteristic of the triplet states. Since it transforms 



IV. EXPERIMENTAL SIGNATURES 

We now turn to the question of the experimental signa- 
tures of such states. Since the order parameter changes 
sign as the Fermi surface is circled, there is no net CDW 
or SDW order which could be measured in, for instance. 



'One can expect, on general grounds, that incommensurate 
singlet or triplet p- or d-wave density-wave order at wavevec- 
tor Q will induce CDW order at 2Q since a term of the form 
^Qp2Q or $Q • $qP2Q is allowed by symmetry in the effective 
action. Nevertheless, we may wish to distinguish such a state 
from one which has only CDW order. 



6 



non-trivially under the point group of the square lattice, 
the triplet order parameter $q will not couple to pho- 
tons, neutrons, or nuclear spins according to ^q-F, where 
F is, respectively, Sn, or /. Said more physically, the 
triplet ordered states do not have anomalous expectation 
values for the spin density but, rather, for spin currents; 
spin currents do not couple simply to these probes. How- 
ever, the order parameter <I>q will couple to such probes 
at second order since its square transforms trivially under 
the point group. Such a coupling will be of the form: 



probe 



(41) 



In the case of photons, this will lead to 2-magnon Raman 
scattering. The coupling to nuclear spins couples directly 
to the nuclear quadrupole moment, and will lead to a 
shift in the nuclear quadrupole resonance frequencies. 

In the presence of disorder, rotational symmetry will 
be broken. Hence, there will be a small coupling, propor- 
tional to the disorder strength, of the Goldstone bosons 
to s-wave probes such as NMR and neutron scattering. 



V. GAPLESS FERMIONIC EXCITATIONS 



The mean-field Hamiltonian is 



H = 



.z. (27r)2 



g<fQfik)4'"\k + Q)4,^{k) 



(42) 



If we define the four component object XAa according to 

(Sl) = (^atfc + g)) ^^^^ 

then the mean-field Hamiltonian can be written in the 
form: 



H 



dPk 
R.B.z. (27r) 



2 X"t(fc) i(e(fc)-e(fc-fQ))T3 -I- 



A(fc)T, + -(e(fc)+e(fc-fg)) (44) 



The integral is over the reduced Brillouin zone. The r's 
are Pauli matrices; the 'flavor' index A= 1,2 on which 
they act has been suppressed. A(fc) is defined by. 



A(fc) = Ao/(fc) = .9<i>Q/(fc) 
The single-quasiparticle energies are: 



E±{k) = -(e(fc) + e(fc + Q)) ± 



(45) 



y{e{k)-t{k + Q)f+A^^{k) (46) 

Let's consider the situation in which there is a node, 
i.e. when the argument of the square root vanishes (we 
discuss below the conditions under which this occurs). 
For simplicity, we will consider the commensurate Q = 
(tt/o, 7r/a) singlet Px density- wave state in a model with 
anisotropic nearest-neighbor hopping: 

e(fc) = — 2t (r cos k^a + cos kyo) (47) 

with r < 1. The mean-field quasiparticle energies are: 

E{k) = ±\J W^{r cos k^a + cofikyof + Ao^sin^fc^ja (48) 

There is a node at kx — 0, kyU — arccos(— r). Expanding 
about this node, 



E{q)^±^vlql+vlql (49) 
with momenta q now measured from the node and 

Vx = Aoa , Vy = 2tay/l — (50) 

The effective Lagrangian for the quasiparticlcs near the 
nodes can be written: 



-Coff = X"^ {dr - T^Vyidy - TxVxidx) Xc 



(51) 



Terms which break the nesting of the Fermi surface, such 
as the chemical potential or next-neighbor hopping, open 
hole pockets at the nodes: 



•^/J — f^X ^Xa 



(52) 



We now turn to the question of when a p- or d-wave 
density wave will have nodal excitations. Let's again be- 
gin with the commensurate Q = (7r/a,7r/a) singlet Px 
density- wave state: 

{^l;''Hk + Q,t)4,pik,t)) ^i\<S>Q\smkxa S"^ (53) 

in a system in which the Fermi surface is nested at Q. 
This state will have gapless excitations if the nodal line 
kx = crosses the Fermi surface. For an open Fermi sur- 
face, this need not be the case. In an anisotropic nearest- 
neighbor tight-binding model, Eq. (^), with r > 1, 
the Fermi surface at half-filling is an open Fermi sur- 
face which does not cross the line kx = 0- Consequently, 
there are no gapless excitations. For r < 1, however, the 
Fermi surface does cross the line kx = 0, and there are 
gapless excitations. 

Are there any thermodynamic singularities at the tran- 
sition at which gapless excitations occur? To answer this 
question, let us consider the mean-field ground state en- 
ergy: 



7 



Eo 



_fk_ 

R.B.Z.(2^ 

(Pk 



R.B.Z 



E{k) 

Ve2(fc) + A2(fc) 



(54) 



The first and second derivatives of Eq are continuous. 
However, tlie tliird derivative of the ground state energy 
with respect to r contains a term of the form: 



(Pk 8e(fc) t^cos^fcajfl 



R.B.Z 



(27r)2 (£2(fc) 



A2(fc)) 



3/2 



(55) 



This term diverges if there is a node on the Fermi surface, 
but is finite otherwise. Hence, the phase with a node on 
the Fermi surface is a critical line with a singular third 
derivative of the ground state energy. We wiU caU this 
phase the 'critical phase' of the Px density wave. Note 
that the second derivative of the ground state energy 
is everywhere continuous but nowhere differentiable in 
the critical phase. It is separated by a third-order phase 
transition from the phase with no gapless excitations, the 
non-critical phase of the Px density wave. 

How does this observation generalize (a) away from 
half-filling and to non-nested Fermi surfaces; and (b) to 
d-wave and/or incommensurate ordering? To answer (a), 
let's change the chemical potential in order to move away 
from half-filling. Now, 



d^Efi _ r (Pk 



Qj,3 



8e(fc) t^cos^kxa 
(27r)2 (e2(fc) + A2(fc))^/^ 



(56) 



Below half-filling, the denominator never diverges. 
Hence, the system is always in the non-critical phase, 
despite the fact that there are gapless excitations. As 
the chemical potential is increased, the system crosses a 
third-order phase transition and enters the critical phase. 
Above half-filling, it is always in the critical phase. 

Suppose, now, that we allow next-nearest neighbor 
hopping t' , thereby spoiling nesting. The ground state 
energy is given by 



En = 



<Pk 
Pk 



E{k) 



E(k)<ti 



(2^)2 



{e{k) + e{k + Q))- 



y{e{k)-e{k + Q)f+AA^{k) 



(57) 



and 



Pk 



4e(fc) t^cos^kxa 



E(fc)<M (2^)' _ ^ Q))2 ^ ^^2(fc) 



3/2 



This diverges if the nodal line of A(fc) crosses the curve 
e{k) = e(fc -|- Q) and this crossing point lies below the 



chemical potential. In such a case, the system is in the 
critical phase. Regardless of the details of the band struc- 
ture, the curve e(fc) — e(fc + Q) is determined by symme- 
try for commensurate Q: it is the set of points for which 
k and k + Q are related by a symmetry of the square 
lattice. For Q = (7r/a,0), e(fc) = e{k + Q) if fc^ ±Tr/2a. 
For Q = (tt/c, 7r/a), e(fc) = e{k + Q) if kx i ky = ±7r/a. 
A d-wave density- wave will always have nodal lines which 
cross the curve e{k) = e(fc + Q)\ a p-wave density wave 
may or may not. If there is no crossing point, or the cross- 
ing point is not below the chemical potential, the system 
is in the non-critical phase. Again, the non-critical phase 
can have gapless excitations. 

In the case of incommensurate ordering, similar con- 
siderations hold. Let us suppose that the Fermi surface 
is nested at incommensurate Q, i.e. if k is on the Fermi 
surface, then fc -I- Q is as well, and e(fc) = e{k + Q) — fi. 
If the Fermi surface intersects the nodal lines of A(fc), 
then there will be gapless nodal excitations. If the chem- 
ical potential is now lowered or the hopping parameters 
are changed, so that the Fermi surface is no longer per- 
fectly nested, then the nodes will open into hole pockets. 
Again, as the nesting condition is approached, a third- 
order phase transition will occur, as in the commensurate 
case. 

In summary, the system will be in a 'critical' state if 
nodal points are at or below the Fermi surface. Oth- 
erwise, the system will be 'non-critical', whether or not 
there other gapless excitations. The transition between 
these two states and the entire critical phase is character- 
ized, in mean-field-theory, by a diverent third derivative 
of the ground state energy. There is no reason to mistrust 
mean-field theory since there aren't strong order parame- 
ter fluctuations which might destabilize out calculations. 



VI. TRANSITIONS TO SUPERCONDUCTING 
STATES 

As Zhang |^ has recently emphasized, enhanced sym- 
metry can be dynamically generated at a critical point 
between two different ordered electronic states. The fo- 
cus of that work was a critical point between an antifer- 
romagnet and a dx2-.y2 superconductor. In earlier work, 
Yang identified an SU{2) symmetry (which, together 
with SU{2) spin-rotational symmetry, trivially forms an 
0(4) = SU{2) X SU{2) X Z2) which is an exact symmetry 
of the Hubbard model at half-filling with /i = U/2. This 
symmetry would be dynamically generated at a critical 
point between a CDW and an s-wave superconductor. 
We now consider the modification of this idea to p- and 
d-wave ordering. 

We first consider a transition at half-filling between a 
singlet commensurate dx2_y2 density- wave and a dx^-yi. 
superconductor. We group the two order parameters into 
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a vector 



f{k) = 



( N/2Re|('^|(fc+§)v4(- 
V2Im{^|(fc + f)^j(- 



V i(V"t(A; + Q+f)V^„(fc- §)) / 




(58) 



if Hint is given by 
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,(0,1) x(04)^ 



If, following Yang g, we introduce the following SU{2) 
generators which we will call pseudospin SU{2) 



(1,0), 



,(l,0);(i.O)(^),(l,0)(^) 



Al"-^(q)AlV(q)- 



Q (9)- 

(0,1) 

Q 



0+ = 



R.B.Z, 



(27r)2 



R.B.Z. (27r)2 



R.B.Z. 



(27r)2 



V'"t(fc + Q)Va(fc + g) 

zV|(fc)^|(-fc + Q) 
«V'T(fc)V'i(-A: + g) 



where 

A(o.") = 



A, 



(1,0) 



(59) 



then the order parameters form a triplet under this 
SU{2), 



$+(<?) /(A) \ 
*o(9)/(fc) 
/(fc)/ 



;^|(fc+i)v^I(-fc+§) 

i(V"t(fc + Q+f)'^a(fc- § 

(^T(fc + §)V^i(-fc + §)) 



(60) 



>(i,o) 



U0,1) ^ 



(2^ 
(2^ 



(2^ 



(2^ 



/(fc)* 



1*? 



/(fc)e"'=**ca (fc+l) 



(-1) 



S/3 



^Ap (k - 



There is a small but important difference between our 
pseudospin SU{2) and Yang's the factors of i in the 
definitions of O^. These are necessary since a commen- 
surate d^2_y2 density- wave breaks T, while a dj.2_y2 su- 
perconductor does not. Consequently, our pseudospin 
SU{2) does not commute with T, which is an inversion 
followed by a rotation by tt about the 3-axis. 

The electron fields transform as a doublet under the 
pseudospin SU{2) as well as the spin SU{2). We will 
group them into 4-component objects "^Aa, where A is 
the pseudospin index, A = 1,2, and a is the spin index. 



.(0,1) 



i(i.i) 



i(i^i) 



_fk_ 
(2^ 



_fk_ 



(2^ 



■/(A)e-^^vl/ 



4 

" 2 

./3 



fc- 



/(fc)* 



C7 



B/3 



(64) 



*la 
*2a 



■00 (fc) 



Q) 



(61) 



A 'microscopic' Hamiltonian which is 0(4) invariant 
can be written down: 



H = Ho 
d^k 



Hi, 



R.B 



.z. (2^)= 



^{k)^ 



Aaf 



Act 



(62) 
(63) 



These 'microscopic' Hamiltonians describe electrons 
at half-filling with a nested Fermi surface and interac- 
tions which favor density-wave and superconducting or- 
der equally. In other words, they describe a critical point 
at half- filling between a dx2-y2 density- wave and a dx2-y2 
superconductor. Near the critical point, we can focus on 
the low-energy degrees of freedom: the Goldstone modes 
and the nodal fermionic excitations. We can write down 
an 0(4) invariant action for this: 



^cff = / dr 



^We will use underlined lowercase Roman letters such as 
i = 1, 2, 3 to denote pseudospin triplet indices and uppercase 
Roman letters to denote peudospin doublet indices A = 1,2. 
Lowercase Roman indices a = 1,2,3 will be vector indices 
(i.e. real spin triplet indices) and Greek letters a = 1,2 will 
be used for real spin SU{2) spinor indices. Pauli matrices r- 
will be used for pseudospin, while cr" will be reserved for spin. 



«5 



dr ■ 



d^fc 

d^k 



Aa^ 



d^q 



(27r)2 (27r)2 



(fc) {dr-e{k))-^Aa{k) 
Hi) f(k) X 



**We have only written down the quartic terms; higher-order 
0(4) invariants also exist, but they are irrelevant at weak 
coupling. 
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ea;3*^"^ (k 



iBBC 



1 1 2 



$1 



$2 _ 



(66) 



(67) 



transitions 



spin-triplets. Hence, the effective field theory for such a 
transition takes the form: 



ScS — 

(65) 



In this Lagrangian, we have rescaled all of the veloci- 
ties and stiffnesses to 1. In general, these quantities will 
be different - breaking the 0(4) symmetry - and this 
cannot be done. Symmety-breaking terms will be briefly 
addressed below. 

The transition between the d^^-y^ density- wave and 
the dx^-y2 superconductor is driven by a pseudospin-2 
symmetry-breaking field, which we will call u. 



For u < 0, the 3-axis is an easy axis and the dx2^y2 
density-wave state is favored; for u > 0, the 1 — 2-plane 
is an easy plane and the dx2-y2 superconducting state is 
favored. 

We can move away from a nested Fermi surface by tun- 
ing the chemical potential or a next-neighbor hopping pa- 
rameter. Such effects are encapsulated by a pseudospin-1 
symmetry-breaking term: 



where O'^ is the pseudospin SU{2) generator defined 
above. If u = 0, will immediately force the pseudospin 
into the 1 — 2 plane - i.e. the supercondcutor will be fa- 
vored. If 7i < 0, the dx2-y2 density-wave state will be fa- 
vored until fic oc {\/—u). At this point, a first-order phase 
transition - the pseudospin-flop transition - will occur at 
which the pseudospin switches from an easy-axis phase 
to an easy-plane phase. If we allow <I>o to have a different 
velocity than $±, then this first order phase transition 
can become two second order phase transitions. Depend- 
ing on the values of these parameters and the strength of 
quantum fluctuations, the intervening phase can either 
have both types of order or neither. 

The critical point occurs when the jump in $ is tuned 
to zero. Hence, it is a tricritical point. At such a critical 
point, 0(4)-breaking terms can scale to zero. The criti- 
cal point and the quantum critical region j2j,|2^ are de- 
scribed by the physics of the critical fluctuations coupled 
to nodal fermionic excitations. By arguments similar to 
those of |^l|, the nodal fermions are neutral, spin-1/2 ob- 
jects. A more detailed analysis will be given elsewhere 

Similar conclusions can be drawn for dr^y and d^2_y2 + 



the latter case is particularly simple 
since there are no fermions. In the case of transi- 
tions between p3;-wave density-wave and superconduct- 
ing states, the order parameters are both pseudospin and 



dT 



d^k 
{2ny 



(a,-e(fc))*Aa + 



«5 



a^Jep^^^'^^ [-k + I) rl^e^'^M/^^t 



-k + ^ 



drd^x 



(a^%)V-r<i>,?a>^ 



4! 



where 



$1 

$3 



^ V2Rc|/v4(fc)e^"<^^;^(-fc)U \ 

V2Im |(Jv:i(fc) e''"<'^ V^(-fc) , 
V i{^k^\k + Q)alP^p{k)) 



(68) 



VII. DISCUSSION 

In this paper, we have discussed the properties of or- 
dered states in which particle-hole pairs with non-zero 
angular momentum condense. These states generalize 
charge- and spin-density wave states in the same way 
that p- and d-wave superconductors generalize s-wave su- 
perconductivity. However, unlike in the superconducting 
case - where the Meissner effect follows directly from 
the broken symmetry, irrespective of the pairing channel 
- the angular variation of the condensate makes p- and 
d-wave density- wave ordering difhcult to detect. Exper- 
iments seeking to uncover such order must (a) be sen- 
sitive to spatial variations of kinetic energy or currents 
or (b) measure higher-order correlations of the charge 
or spin density. We explained how /iSR, neutron scat- 
tering, NQR, and Raman scattering can be used in this 
regard. Impurities, which break rotational invariance, 
would cause admixture of p- or d-wave ordering with s- 
wave ordering. It is natural to wonder whether experi- 
ments which appear to detect SDW order should be re- 
examined to see if they have actually uncovered p- or 
d-wave order which, as a result of impurities, is mas- 
querading as s-wave order. 

As in the superconducting case, the non-trivial pairing 
symmetry can lead to the existence of nodal excitations. 
As parameters such as the chemical potential or next- 
neighbor hopping are varied, nodal excitations appear at 
a transition which is third-order in mean field theory. 
The 'phase' with nodal excitations is always critical. 

The analogies between p- and d-wave density- wave or- 
dering and p- and d-wave superconductivity begs the 
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question: what is the nature of a phase transition be- 
tween such states? In answering this question, we are 
led to one of the motivations of this work. The pseudo- 
gap regime of the cuprate superconductors exhibits some 
properties which can be associated with d^i-yi order- 
ing. One explanation is that some features of the (1^2 _y2 
superconducting state have been inherited. However, it 
is natural to inquire whether the physics of this regime 
could also be determined in part by proximity to a dx2-y2 
density-wave state or the transition between the density- 
wave and superconducting states. In other words, we ask 
whether the physics of the pseudo-gap regime should be 
described by a theory which incorporates fluctuations be- 
tween dx2-y2 density- wave and superconducting states. 
In this connection, we note that the physics of the critical 
point between dr^2_y2 density- wave and superconducting 
states bears a rough resemblance to that of the SU{2) 
mean field theory of the t~ J model . In that the- 

ory, the gauge field parametrizes fluctuations between the 
dx2-y2 density- wave and superconducting states, a role 
played in our analysis by the Goldstone bosons of the 
0(4) effective theory. We also note that the Nodal Liquid 
state pl|-p3t shares many features of the dx2_y2 density- 
wave and superconducting states. These issues and their 
possible relevance to the cuprates will be further explored 
elsewhere. 

I would like to thank Sudip Chakravarty, Stuart 
Brown, Subir Sachdev, and, especially, Dror Orgad for 
discussions. 
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